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Abstract 

We study the dynamics of a spatially inhomogeneous quantum Xcf)^ field 
theory in 1 + 1 dimensions in the Hartree approximation. In particular, we 
investigate the long-time behavior of this approximation in a variety of con- 
trolled situations, both at zero and finite temperature. The observed behavior 
is much richer than that in the spatially homogeneous case. Nevertheless, we 
show that the fields fail to thermalize in a canonical sense, as expected from 
analogous results in closely related (mean field) transport theory. We ar- 
gue that this dynamical approximation is best suited as a means to study 
the short-time decay of spatially inhomogeneous fields and in the dynamics 
of coherent quasi-classical inhomogeneous configurations (e.g. solitons) in a 
background of dynamical self-consistent quantum fluctuations. 



PACS Numbers: 03.70. +k, 05.70.Ln, ll.lO.-z, 11.15.Kc 
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I. INTRODUCTION 



In recent years much progress has been made in understanding the real time dynamics of 
spatially homogeneous quantum field theories. This has been made possible by considering 
mean field approximations {e.g., Hartree or the large- limit) to the full quantum dynamics 
of these systems. The advantage of such techniques is that they permit the approximate 
evolution of an arbitrarily out-of-equilibrium quantum field theory under many different 
circumstances, such as in the presence of large mean fields and/or at symmetry breaking 
phase transitions, where perturbative methods fail. 

Numerous applications of this type of dynamics have been considered, ranging from 
symmetry breaking phase transitions decoherence and dephasing 0, the formation of 
disoriented chiral condensates at the collisions of relativistic heavy ions to the re-heating 
and preheating of the early universe [Q, to name just a few. Nevertheless, the approximations 
involved in this kind of evolution lead to two serious limitations, standing in the way of 
making the results apply to more general and realistic circumstances. The first drawback 
is that the long-time behavior of the system is incorrect due to the absence of collisional 
effects. The second is that only spatially homogeneous situations could, until recently [^Q, 
be considered. Both limitations are not a matter of principle; they can be overcome with 
large increases in computational effort, although theoretical progress, e.g., the development 
of a time-local formulation of the quantum field dynamics including scattering, would be 
welcomed. 

In this paper we lift the constraint on the spatial homogeneity of the system. Mean field 
theories can be very naturally generalized to spatially inhomogeneous circumstances, allow- 
ing the study of the (approximate) evolution of quantum fiuctuations in the background of 
dynamical, spatially dependent fields. This generalization results in a demanding numerical 
problem, requiring the solution of a large set of coupled partial differential equations. As 
will be made clear later, the size of the computational effort scales as L^^, where L is the 
linear dimension of the problem and D the number of space dimensions. Thus problems in 
lower spatial dimensions are significantly more accessible. 

One question of great general interest is whether by considering spatially inhomogeneous 
situations one improves on the long-time properties of the mean field approximation. In 
particular we will investigate if the long-time regime is in any sense universal and, in par- 
ticular, if it corresponds to canonical thermal equilibrium. Salle, Smit and Vink ^ have 
recently found that, under averaging over an ensemble of certain classes of initial conditions 
for the mean field, the statistical occupation number of the fiuctuations may dynamically 
approximate a Bose-Einstein distribution over a limited range of times. Although these re- 
sults are very interesting, we believe that the averaging over ensembles of initial mean fields 
is not justified from first principles, since the mean field is by definition already the quan- 
tum ensemble average one-point function. The potential power of inhomogeneous quantum 
evolutions is to handle situations where a specific spatially dependent expectation value of 
the quantum field exists due to some perturbation of thermal equilibrium or of another 
asymptoticly stable state. Restrictions to particular types of initial conditions defeat the 
general purpose of the method. For these reasons we will adopt a more general stance here 
and investigate the main characteristics of the mean field dynamics without any additional 
ensemble averaging. 
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We do this in the context of a relativistic quantum field theory in 1 + 1 dimensions. 
In section |l| we define the quantum field theory, and the mean field equations of motion. In 
section ^ we make our choice of initial conditions and discuss the renormalization scheme 
necessary to make the theory independent of the ultraviolet cutoff. We proceed in section |V 
to analyze the decay of a family of Gaussian spatial profiles. This allows us to investigate the 
dynamical response of the quantum fiuctuation to this field background. We characterize 
general properties of the dynamics of the mean field and of the quantum fiuctuations and 
pay particular attention to their long-time limit. We then proceed in section 0, in the spirit 
of transport theory, to study the evolution of a mean field small perturbation to a state of 
thermal equilibrium. Together the results of section |rV| and [V] allow us to determine whether 
the model thermalizes at long times in circumstances of general interest. Finally we present 
our conclusion and outlook for further applications. 



II. THEORY AND DEFINITIONS 

We study the dynamics of a relativistic quantum theory in 1 + 1 dimensions. We 
write the Lagrangian density as 

1 Tfi ^ \ 

^ = -^{d,m'ct^)-^^''\^'- (1) 

The field is taken to be real. We will work in units in which the speed of light c and 
Boltzmann's constant are unity c = 1, ks = ^■ 

There are several different ways of obtaining mean field approximations to the quantum 
field theory, by writing truncated effective actions (e.g. in the large limit [0])! by making 
a Gaussian variational ansatz in the context of a Schrodinger functional approach ||^, or 
by explicitly expressing higher point field correlation functions in terms of the one- and 
two-point functions, thus effectively truncating the Dyson-Schwinger hierarchy. While these 
procedures do not lead to the same approximation to the full thermo(dynamics), the resulting 
equations share much in common. Below we adopt the familiar Hartree approximation to 
the dynamical equations for the one- and two-point functions of the quantum field (f). The 
type of evolution described in this paper generalizes to any other mean field approximation 
with a similar amount of computational effort. 

Following these remarks, we separate the quantum field in a mean field (/9 = (0) and a 
fiuctuation field ip {{ip) = 0), such that 

^{x,t) = ip{x,t) +ij{x,t). (2) 

Our mean field approximation (which coincides with the familiar Hartree approximation) 
corresponds to keeping only the connected two-point function in ip. Higher point connected 
correlation functions of the fiuctuations are assumed to be zero. Thus we have, e.g., {tp"^) = 
3(?/'^)^. Equations of motion for the mean field and the Wightman two-point function 

G^{x,t;x',t') = {ip{x,t)ilj{x' ,t')) then become 

□ + + Xif^{x, t) + 3AG>(x, t)] if{x, t) = 0, (3) 
[□ + x{x, t)] G>{x, t- x\ t') = 0, x{x, t) = m^ + 3X<^\x, t) + 3XG>{x, t), 
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with U = d'i-dl, G>{x,t;x',t') = {iP{x,t)iij{x' ,t')) and G>{x,t) = G>{x = x',t = t'). 

< 

The relation between the (bosonic) Wightman functions and the time ordered (Feyn- 
man) two-point function Gpixjt] x' ,t') is 

G>{x,t;x',t') = G<{x',t';x,t), G>{x,t;x' ,t') = G<*{x,t;x' ,t'), (4) 
Gf{x, t; x', t') = 0(t - t')G>{x, t; x', f) + Q{t' - t)G<{x, t; x', t'), 

where 6(t) is the usual step function. The advantage of working with the Wightman func- 

< 

tions G^ is that they obey homogeneous equations of motion. 

It is convenient to solve Eqs. (^) by decomposing G^ {x,t; x' ,t') in a complete orthonor- 
mal basis of mode fields, which we shall denote by {ijjk{x,t)}. The mode index k may be 
continuous or discrete. In a general context we will denote traces over the basis of field 
modes as sums over k. Then we can write the quantum field ip{x, t) as an expansion in this 
basis: 

^(a^) t) = Y^ [aki^kix, t) + ali)l{x, t)} , (5) 
k 

where a^, a\ are the annihilation and creation operators, respectively, obeying canonical 
commutation relations 



hq, K, CLq] = al, al = 0. (6) 



Using (H) and (D and {a\aq) = dtqUsik) we can write G^ {x,t; x' ,t') as 

G>{x, t- x\ t') = {{uBik) + t') + nB{k)ijl{x, t)M^', 0} , 

k 

where nsik) is the Bose-Einstein thermal occupation number distribution 

with ujk the energy eigenvalue associated with the eigenfield ipk{x,t). In particular, the 
equal-point function G-^ {x, t) acquires the familiar form 

G>{x,t) = J2 l^Mk) + l]M^,t)rk{x,t). (8) 

k 

The spectrum reduces to the vacuum as T — > 0, since ^ 0. 

The equation of motion for G^{x,t;x',t') can now be written in terms of the mode 
functions ipk- 

p + x{x,t)]Mx,t) = 0, (9) 

which must be solved for all k. Additionally x{x, t) can be expressed entirely in terms of ip 
and tpk through Eq. (|^), so that we obtain a set of self-consistent partial differential equations 
for the fields ip{x,t) and for each of the ipkixit). This set is formally infinite; to obtain a 
good approximation to the dynamics in the continuum we must take a very large number 
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of mode fields ilJk{x,t), thus the problem becomes demanding in terms of computational 
resources. 

In practice, the fields (p{x,t),ipk{x,t) are evolved on a spatial lattice. Renormalizability 
requires that the mode fields ipk must reduce to the vacuum at large k and become plane 
waves. This limit places a restriction on the modes that can be resolved dynamically: 
for a spatial lattice of sites, with spacing Ax, modes with wavenumber k > 27c/Ax 
cannot be kept as dynamical fields. Using periodic boundary conditions in space, the nth 
wavenumber is fc„ = 2im/N /\x. Thus the number of mode functions is restricted to 
per linear dimension. Then we see that the computational effort in D space dimensions 
corresponds to that of solving A^^ coupled partial differential equations on a grid of size A^^ 
points, i.e. it scales as A^^^ per time step, or like the D-dimensional volume squared. This 
scaling makes the evolution on large spatial lattices in three dimensions very demanding^. 
For the rest of this paper we will be concerned with time evolution of one-dimensional fields. 

Finally, we need to specify a basis {^^(a^, ^)}, together with the scalar product 

t).^pq{x, t)) = z y dx {^lix, t)dt^lJq{x, t) - iljq{x, t)dtipl{x, t)} . (10) 

It follows directly from the form of the equations of motion (|^) and that the scalar 
product ([10|) is invariant under the dynamics, i.e., dt {ipk{x,t).xl'q{x,t)) = for all k,q. This 
guarantees that the orthonormality of the basis, chosen at the initial time, remains valid 
to all subsequent times, although the set of fields ipk will evolve away from their initial 
conditions. 

At the initial time, it is convenient to parameterize the set of functions ipk by 



Mx,t) = ^^^e-^-^'g,{x). (11) 

The canonical equal-time commutation relations for the fundamental quantum field '?/^(x, t) 

( [V^(x, t), df^ix', t)] ) = ih6^{x - x') (12) 
must hold. This places a further requirement that the Wronskian condition 

ijjkix, t)dt'il>l{x', i) - dti>kix, i)iljl{x\ t) = ihgk{x)gl{x'), (13) 
together with the completeness condition 

Y.9k{x)gl{x') = 5''{x-x'), (14) 

k 

must be satisfied. In the spatially homogeneous case one chooses gk{x) to be plane waves 
gk{x) = e*^^, which clearly satisfy (p!3| ) and (|14|). Once satisfied at t = 0, the stationarity of 



the scalar product also guarantees these properties at later times. 



The limits of present computational power allow for the solution of the three dimensional problem 
on a lattice of size N^, with ~ O(IO^), for short times. 
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Given ([TOD-(p!lD, the purely spatial fields gk{x) satisfy the static eigenvalue problem 



-V^. + x{x)\ gk{x) = uj^gkix). (15) 
The scalar product (0) can now be written in terms of Qk- 

i9kix).gg{x)) = dx gl{x)gg{x), (16) 



which is familiar from quantum mechanics. With these definitions orthogonality of the set 
{gk} is equivalent to orthogonality in {ipk}- 

III. INITIAL CONDITIONS AND RENORMALIZATION 

To fully specify our initial value problem we need to choose initial conditions for the 
fields. The initial fields must also be spatially periodic. For the mean field we choose a 
family of Gaussian profiles 



ip{x, t = 0) = (pQ exp 



'2A 



(17) 



where (po, A are real parameters. In addition we choose the conjugate momentum dtip{x, t = 
0) = 0, for simplicity. The choice of the set ipk carries some arbitrarily for small k. As we 
have seen above, only a self-consistent set of modes solved in the background of (|l3) can 
supply us with a globally static solution. This type of solution is currently being investigated 
jOIJTOd for specific profiles of ip. 

The simplest possible set of mode fields satisfying all necessary requirements are plane 
waves 

gk{x) = e'"-'. (18) 
The periodicity of the spatial boundary conditions enforces k ^ kn = 27m/ L. For numerical 



work we take the corresponding lattice solution with c<j„ = ym^ + 2(1 — cos(A;„Ax))/Ax^, 
where L = NAx is the physical size of the system, N its number of sites and Ax the lattice 
spacing. 

Like most quantum field theories, our model as it stands is unphysical, as it collects 
infinite vacuum contributions as Ax — >■ 0. Simple power counting shows that in 1 + 1 
dimensions is superrenormalizable: there is a single 1-loop logarithmic ultraviolet diver- 
gence in the self-energy. There is also an additional quadratic divergence in the energy, 
characteristic of the free theory. 

To see this, consider the vacuum solutions to equation (|]). For the purposes of renor- 
malization we can assume ip to be spatially homogeneous. We write x 

3A 

X = ml + 3Xip^ + — dk \i;km2nBik) + I] , (19) 

TC Jo 

where the bare mass is understood to be cutoff dependent in a way that will cancel the 
divergence resulting from the integral. 
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To define the spectrum of vacuum excitations, we clioose x to be tlie mass squared 
of the excitations, obtained from the curvature of the classical interaction potential at its 
m? . Then, without symmetry breaking y9 = 0, we obtain the renormalization 



mmimum, x 
condition 



m 



m 



+ — I dk 



^ ' TT Jo 2Uk^ 



with iUk = i/fc^ + X = + m'^- This results in 



ml 



X 



vr Jo 2^/k'^ + m?' 

3A { 

+ 3A(^2 + — / dkl \^k? \^nB{k) 

TT JO [ 



(20) 



(21) 
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where the last relation is valid at all times. Clearly we need to take A <C T to obtain results 
that are independent of A. 

The choice of plane waves for {(^^(a;)}, apart from its simplicity, is not ideal, in the 
sense that the introduction of an inhomogeneous mean field causes the initial choice of 
basis to shift quickly to adapt to this background. This leads to a fast transient in the 
spatial profile of the mode functions ipk over the set of long length scales that characterize 
the spatial inhomogeneity of the mean field. Consequently, this effect does not change the 
renormalization, which is an ultraviolet property. 

When computing the energy, we have, in addition to the logarithmic self-energy diver- 
gence, a quadratic divergence characteristic of the free theory. The total energy has the 
form 



E 



dx {\dt(p{x)\^ + \V(p{x)\^ + m^ip{x) 



+ [m' + 3Xcp' + ^AG>(x)) G>ix) + [|V^,(x)p + \dtMx)\' 



(22) 



The quadratic divergence arises from the last integral in 
subtract the vacuum contribution to this term 



AE 



7r Jo 



dk 



k^ + ujl 
2uJk 



7r 



dk k 



To remove it we need to 



(23) 



which is clearly a quadratically divergent integral with the ultraviolet cutoff A. To render 
the energy finite it is sufficient to subtract the ultraviolet behavior of (^). By subtracting 
the full integral in (^3]) we impose the renormalization condition that the physical energy at 
zero temperature is contained exclusively in the profile of the mean field (f. 

We are now ready to study the non-equilibrium evolution of the quantum fields. 



IV. THE TIME EVOLUTION OF INHOMOGENEOUS MEAN FIELDS 

In this section we discuss the dynamical properties of the spatially inhomogeneous 
Hartree approximation. We solve Eqs. (H) and (^, in the region [— L/2,L/2], with pe- 
riodic boundary condition in space. We discretize the computational domain into + 1 
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time = 50 



time = 50 




points, with xj = jAx, j = -N/2, -N/2 + 1, 0, 1, 2, N/2, where Ax = L/N. We ap- 
proximate G^{x,t) using N complex base fields. A fourth-order symplectic integrator with 
a fixed time step At was used to advance the solution to the final time tf. In the results 
shown below, L = 128, N = 1024, At = 0.0025. 

As a zeroth-order check, we verified that the total energy is conserved to within 0.001%, 
or one part per 100,000. We have also made a "convergence study", where, for one set of 
initial parameters, we performed the simulations using three different spatial discretization 
Ax, and found that the integral features of the solutions are unchanged with refinements of 
spatial resolutions. 

In Figure |l| we show several snapshots of the evolution of an initial Gaussian profile with 
00 = 2.0, A = 1, and with parameters A = 0.1, h = 1, m? = 1. This situation corresponds 
to a macroscopic local over-density of the field 0, which could be achieved by turning on an 
external potential. At t = the external perturbation is released and the over-density is free 
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FIG. 2. Snapshots of the power spectrum of the mean field ip. The power stays approxi- 
mately contained within the initial Gaussian envelope (dot-dashed), while decaying slowly as time 
progresses. 

to propagate and decay. We see that the Gaussian profile oscillates and decays by emitting 
wave packets that travel outwards from the perturbation point, both in its own profile, as it 
would happen in the purely classical case, and in the fluctuation modes, which correspond 
to particles created by the perturbation. 

Because the fluctuations are massive, the medium is dispersive, and different wavelengths 
propagate outwards at different velocities. We observed that their propagation velocity is 
bounded by c = 1, and that shorter wavelengths travel faster, as is apparent in the first 
few frames of Fig. |1|. Meanwhile, some of the energy of the mean field is transfered to the 
fluctuations, which qualitatively display a similar behavior of wave packets traveling outward 
from the perturbation point. For late times the wave packets collide due to the periodic 
nature of the spatial boundary conditions, and continue to propagate back and forth in the 
computational domain. 

An important feature of the evolution is that the mean field does not acquire large power 
over small spatial scales. In fact, as shown in Fig. most of the power spectrum of the mean 
field stays contained within the initial Gaussian envelope, although there is some observable 
enhancement over smaller wavelengths due to the back reaction of the fluctuations. 

We now turn to the question of the long-time behavior of the dynamics. As a preamble 
we should keep in mind that our approximation resumes interaction effects that appear 
at 1-loop order ~ A. Consequently we can expect the approximation to be valid up to 
times t ~ 1/A = 10, for the case studies shown hereQ. CoUisional effects, not included in 



^We performed a large variety of other evolutions, obtained by varying A and m and the initial 
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the Hartree approximation, first appear at two-loops (~ A^) and are necessary to render 
the evolution physically valid to times of order t ~ Thus pushing the mean field 

approximation to large times requires a substantial optimism. Any conclusions it suggests 
for the properties of the physical system it attempts to describe in this regime should be 
met with a healthy dose of skepticism. 

In any case, as we have discussed above, as time progresses the mean field loses some of 
its initial energy to the fluctuations. These, however, do not always carry it away efficiently 
(see also next section), leading at late times to a local dynamical balance between mean 
field and fluctuations that does not result in the total decay of the former. In this sense we 
see already that a state of canonical equilibrium is not reached, since the mean field remains 
non-vanishing. Temporal averages of local mean fields, however, will yield approximately 
zero. In this sense we may think that the system has reached a state of microcanonical 
equilibrium and that any local observable, averaged over time, may mimic a corresponding 
average obtained from the canonical ensemble. 

To help us decide whether we may have reached a state of microcanonical equilibrium, we 
analyze the spectrum of fluctuations ipk- Again, a canonical system in contact with a thermal 
bath would be expected to thermalize (in the mean field approximation) at a temperature 
consistent with equipartition. If the system thermalizes, the fluctuation spectrum would 
satisfy 

im') = ^ [2nB{uJu) + 1] , (iTTfeH = ^ [2nB{iVk) + 1] , (24) 



where ujk = \J k"^ + m?{T), or its form on the spatial lattice, and 7rfe(a;, t) = dtipki^, t)- These 
reduce to their classical Boltzmann forms in the limit of — 

{M = ^, (|vr.P)=T. (25) 

Because of the exponential decay of nB{uJk) for large uj^, the vacuum contribution is typically 
dominant with T = initial conditions. We plot instead {\ipk\) — {\'<Pk\)o, by subtracting out 
the vacuum piece. 

To investigate whether the late time fluctuation power spectrum is in any sense universal 
we set up a family of initial conditions by varying 0o and A under the restriction that the 
renormalized total energy of the system E stays constant, keeping A = 0.1 fixed. Table I 
shows the values of 0o aiid A used, for two different values of the energy E. Fig. |^ shows 
several late-time fluctuation power spectra for initial Gaussian mean field profiles, while 
Fig. 1^ displays the corresponding spectra of the fluctuation conjugate momenta. After an 
initial transient, the spectra tend to settle to an approximately time independent profile. 
Although some of the general features of the several spectra are similar they differ in, for 
example, amplitude. The fluctuation spectra are reasonably fitted by a power law cu"", with 
2 < a < 4, while that of the conjugate momenta may be fitted by a similar power law, but 
with 1 < a < 2. 



conditions for witli qualitatively similar results to those shown here at A = 0.1 and m = 1. 
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A 


E = 0.008426 


A 


E = 0.855148 


00 


2.00 


0.260128 


2.00 




2.60106 


1.00 


0.200000 


1.00 




2.00000 


0.50 


0.146311 


0.50 




1.46514 


0.25 


0.101603 


0.20 




1.01951 



TABLE I. The height (/)o and variance A of the initial Gaussian mean field (/j(x,0), Eq. 
used to generate a family of evolutions with the same energy E. The results for the late time 
fluctuation power spectrum are shown in Figs. ^ ^. 




FIG. 3. Late-time power spectra of the fluctuations — \il){t = 0)p (with subtracted vacuum) 
for various amplitudes and energies, see Table I. Solid lines show (upper) and uj~^ (lower) 
dependences. 
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E = 0.008426, A = 0.5 



E = 0.008426, A = 2.0 




E = 0.855148, A = 0.5 




E = 0.855148, A = 2.0 




FIG. 4. Late-time power spectra of the fluctuation conjugate momenta l^fV'P — \dtip{t = 0)p 
(with subtracted vacuum contribution) for various amphtudes and energies, see Table I. Sohd hues 
show Lo"^ (upper) and uj~'^ (lower) dependences. 



These results suggest that at late times the system enters into a dynamical steady state 
with fluctuation spectra that are uncharacteristic of canonical thermal equilibrium, either in 
the classical approximation or for quantum fields. Indeed, the characterization of the late- 
time spectrum as an approximate power law with 2 < a < 4 suggests that the spectrum 
of fluctuations is less ultraviolet dominated than in the classical Boltzmann approximation 
{a = 2) but more than in its quantum Bose-Einstein form. We never observed exponential 
suppression of the power spectrum at large u, as would be expected from quantum canonical 
equilibrium. Interestingly, this range of power laws shows that the spectrum of fluctuations 
generated dynamically by the decay of the mean field gives rise to expectation values with 
no divergence in the limit A; — >■ oo. If it holds in higher spatial dimensions this is an 
advantageous feature of the current approximation relative to purely classical field theory. 
Obtaining truly quantum occupation numbers almost certainly requires the inclusion of 



quantum scattering effects |jTT],[T2|. 

The late-time statistical stationarity of the fluctuation power spectra suggests that in its 
long-time limit the fluctuations may decouple from the evolution. This is what happens in 
homogeneous mean field theories, where the long-time limit is characterized by xit) C, 
where C is a constant, dependent on the initial conditions. Then the modes ipk{t) become 
harmonic oscillators with a corrected frequency. The canonical thermal equilibrium behavior 
of x{^} t) clIso demands that it becomes a constant in both space and time. 

In order to probe this behavior, we show the dynamical evolution of xi^^t) ^iid of its 
power spectrum in Fig. ^ At late times the power spectrum of x reaches approximate 
stationarity. Clearly there is very little power for large k, showing that the modes in this 
regime are essentially decoupled. The approximate stationarity of the power spectrum at 
small k suggests, however, that the system has reached a statistical dynamical steady-state 
which does not coincide with canonical thermal equilibrium. 
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FIG. 5. The evolution of xi^j t) for an initial Gaussian mean field profile with (f)Q = 2.0, A = 1, 
A = 0.1, = 1 and k = 1 (left) and of its power spectrum (right). 
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It is clear that initially the mean field induces a strong perturbation in x(a;,t), which 
subsequently decays to smaller values in the central region. For late times xi^y t) remains a 
time-dependent, spatially inhomogeneous field, showing that the system continues to exhibit 
complex dynamics in space-time. This dynamics, however, leads to no significant transport 
of energy among scales, as evidenced by the approximate stationarity of the fiuctuation 
power spectrum. 

Since most non-trivial dynamics happens at early times, we show, in Fig. |^, the very early 
evolution of x{x,t), together with its power spectrum. Recall that a small x corresponds to 
an approximate decoupling of the corresponding mode. As evidenced by Fig. ^ the mean 
field induces initially significant power over a continuum of large scales. This suggests that 
the fluctuations have adapted spatially to the non-trivial profile of (f. The response of the 
high k modes is selective and can be seen to fall into a distribution of bands. These bands 
then split into finer and finer ones to become a continuum at late times. 

A pattern of excited bands in the fiuctuation power spectrum was also observed in 



the case of the spatially homogeneous evolution of similar theories ||T3|. Such bands are 
closely related to exponentially growing solutions of Mathieu's equation, which describes 
the behavior of an oscillator whose frequency is a periodic function, supplied in our case by 
the mean field. In the spatially inhomogeneous case, this pattern is complicated by the fact 
that the mean field possesses a continuum of frequencies each with a different amplitude. 
To sketch these parametric resonances, we neglect the fluctuations in the mode equation 
Eq. (^, since initially this contribution will be small. Then we can write 

a +m'^ + 3Xip'^{x,t)]ipk{x,t) ^0. (26) 

This can be written in Fourier space as 

52 + (^2 ^ ^2)] ^^(^^ t) + 3A I ^^y,(p^)^(p,)V;,(g - - p„ t) = (27) 

The last integral describes mode mixing between the mean fleld and the fluctuation fleld 
ipk- To isolate parametric instabilities we need to further approximate it, by placing pi = 
—p2 = p, which results in 

'd^ + (q^ + m')+3\^']Mq,t)^0, = jdx^\x)= j ^ \v{p)\\ (28) 

This equation can be mapped to Mathieu's equation, 

dt+a + 2e cos t] (q, t) = 0, (29) 

with a = (g2 + + 3A^^(t = 0) /2) /{2mf, e = 3X(p'^{t = 0) /8m^. An instability analysis of 
Mathieu's equation with v^^(t) — v^^(t = O)cos^mt predicts instability bands (exponential 
growing modes) for a = with n an integer. In terms of the Fourier momentum q this 

results in g ~ my/n'^ — 1 = v^, v^, vTS, provided A(^^(t = 0) -C m^. 

Although this analysis of parametric instabilities is necessarily approximate, we observed 
that the values for instability modes predicted by Mathieu's equation are realized to a good 
approximation in the early dynamics either as maxima or minima in the power spectrum 
of fluctuations for k > m. For small fc, however, the population of the power spectrum 
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FIG. 6. Early profiles of xi^T^) (left) and corresponding power spectra (right) for the same 
parameters as in Fig. ^. The higher k modes are populated in bands that later coalesce to form a 
stable continuum spectrum. 
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follows a broad continuum presumably seeded not by parametric instabilities but by terms 
in the mode mixing that act as sources. This behavior is a new feature of the dynamics 
relative to the spatially homogeneous case. Later, more and more bands appear, resulting 
in a complicated behavior, until the spectrum reaches a continuum also for large k. These 
aspects of the evolution are quite complex and lie beyond the scope of our present analysis. 
They will shed light on the detailed mechanisms for energy flow among scales in the model 
at early times, when the mean field approximation to the dynamics of the quantum fields is 
justifiable. 



V. RESPONSE AROUND THERMAL EQUILIBRIUM 

In the previous section we showed some of the general properties of the inhomogeneous 
mean field evolution starting with families of Gaussian profiles for the mean field and the 
quantum fluctuations in vacuum. In this section, in the spirit of transport theory, we 
investigate the properties of the model when perturbed around a self-consistent state of 
thermal equilibrium. 

The mean field equations have as a limiting case a transport mean field theory 
characterized by a Vlasov equation. The derivation of this transport theory is standard 



Tl| and we will not repeat it here. The transport description is appropriate for small soft 
disturbances of the system around a stable asymptotic state, usually thermal equilibrium. 
By solving the corresponding transport theory, one can estimate how the system responds 
over large spatial scales and long times and, in particular, how and if it returns to thermal 
equilibrium. 

In the context of Boltzmann transport theory, the return to thermal equilibrium is re- 
lated to a collision kernel resulting from an imaginary part of the self-energy of the field in 
the medium. For the lowest such self-energy is given by the 2-loop sunset diagram, with 
its imaginary part describing 2^2 particle scattering. Such a self-energy is absent from 
our description. Thus we may simply expect, as is characteristic of mean field transport de- 
scriptions, that disturbances to thermal equilibrium may never decay and that the response 
and final state of the system is described in terms of characteristics of its initial state (such 
as its energy), and is therefore not universal. 

As we have seen in Eq. (P), thermal initial conditions are easily included through the 
definition of G^{x,t). To ensure that we start with a self-consistent (static) solution, we 
must solve the eigenvalue problem corresponding to the state of thermal equilibrium in the 
mean field approximation. In thermal equilibrium cp = and the tadpole diagram acquire a 
momentum-independent mass correction Am^(T). The self-consistent solution for the modes 



i/^k then is in the form of plane waves (|TTD, (|18|) , with frequency uik = yk'^ + m? + Am2(T). 
The mass correction A'm'^{T) obeys the gap equation 



ul = A;2 + m2 + Am^(T), 

2nB{uJk) + 1 



A 2 /-A 

Am^ = / dk 

n Jo 



2ujk 2VA;2 + w? 



(30) 



where we have used a renormalization prescription such that the vacuum tadpole is fully 
subtracted so that the result coincides with (^O]), (|2T|) as T ^ 0. We implement this solution 
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T Am^{T) 

"OO 0.0000 

0.5 0.0116 

1.0 0.0523 

2.0 0.1535 

3.0 0.2577 

5.0 0.4580 

10.0 0.9042 



TABLE II. The thermal mass correction Am^(T) vs. the temperature T. The computations 
used a lattice propagator, and A = 0.1, m = 1. 



on our spatial lattice modes ipk in order to produce a static solution for the numerical 
evolution. Some sample results are shown in Table II. We verified that the solution resulting 
from this procedure is static in the absence of the mean field or an external perturbation. 

An important point to note is that the form of the Bose-Einstein distribution is not 
essential for the stationarity of the solution. Equation (^) can be satisfied for any other 
occupation number, provided the integral over k exists. This contrasts once again with the 
integrals in Boltzmann coUisional kernels, for which the Bose-Einstein distribution is usually 
the unique solution. 

We can now explore the response of the fiuctuations to a self-consistently evolving mean 
field or to a static external perturbation. In the examples below we will work at T = 3, 
which is a high temperature situation since T > m(T) = 1.12. 

Consider first an external field ipext- We assume that a static solution will be reached 
in this external field and parameterize the linear deviation from thermal equilibrium by 
generalizing the equilibrium number distribution to a general slowly varying function of 
space and time n{x, k) such that 

n{x, k) = UBiuJk) + 6{x) ^^^^^''\ (31) 
By requiring stationarity of the solution we can conclude that the field 6{x) is 

,2 



V^ext 

r — 



S{x) = I J , , . (32) 



We can now compare this simple expectation with the result from the evolution. An example 
is shown in Fig.0 for a Gaussian profile of the form (plTl) with 0o = 0.2 and A = \ 

It is striking to see that the fiuctuations do not completely cancel out the externally 
imposed field, and that a globally static solution is not reached at long times. Instead, the 
fiuctuations adapt to the external potential to form a local bound state and some of the 
energy perturbation introduced at t = 0, when the mean field was switched on, is carried 
away by wave packets propagating outwards. 

For the same mean field self-consistently evolving with its fiuctuations the result is rem- 



iniscent of that of section IV, see Fig.H. As in the case of the external perturbation, the 



self-consistently evolving mean field is partially shielded by fiuctuations, but not enough to 
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FIG. 7. G-^{x, t) — (red) showing the response of the fluctuations to an externally imposed 
Gaussian field cp^^^ (blue). Eq. (p!7[), with A = 1 and 0o = 0.2 and at temperature T = 3. 
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FIG. 8. G^{x,t) — (red) showing the response of the fluctuations to a self consistently 
evolving mean field squared (blue) started from a Gaussian profile Eq. (p!7|), with A = 1 

and (/>o = 0.2 and at temperature T = 3 



19 



yield a globally static solution. Instead, as the mean field evolves and its profile changes in 
space-time, the fluctuations track that change on essentially the same time scale after an 
initial fast transient. 

The properties of self-consistently determined fluctuations in fixed external backgrounds 
both at zero and finite temperature can be naturally studied in our dynamical approximation. 
A procedure to compute minimal energy fluctuation modes in these potentials can be easily 
generated by introducing artificial dissipation in the evolution. 



VI. CONCLUSIONS 

We have analyzed the dynamical properties of an inhomogeneous relativistic quantum 
Ac/)'* field theory in 1+1 dimension in the Hartree approximation. We considered the evolution 
of the fields starting both from a vacuum state and from a state of thermal equilibrium at 
finite temperature, perturbed by a soft mean field profile. 

The resulting evolution can be very complex, opening up many new possibilities in the 
dynamics of quantum fields arbitrarily away of equilibrium. Spatially inhomogeneous fields 
are necessary, e.g., in the study of the plasma formed at the ultra-relativistic collision of 
heavy ions and its subsequent anisotropic cooling and in considering the evolution of charged 
quantum fields in a mean magnetic field. 

The transport of energy from the mean field to the fluctuations proceeds at early times 
selectively into intervals of spatial scales, which later coalesce to create a continuous spec- 
trum, which becomes almost static at late times. This process, similarly to the spatially 
homogeneous case, shows effective time irreversibility emerging within the context of a fully 
unitary quantum field evolution. Although some features of the late fluctuation spectrum 
may be quite general, others, such as its amplitude, are dependent on the initial state. The 
dependence of the power spectrum on uj suggests that the mean field evolution generates 
occupation numbers that do not reproduce the quantum behavior of the exact quantum field 
theory but that, at the same time, are better behaved than those of the classical thermal 
theory. In particular this spectrum leads to no new ultraviolet divergence. It will be inter- 
esting if this characteristic holds in higher spatial dimensions where the classical Boltzmann 
distribution leads to ultraviolet divergent thermal contributions to expectation values of field 
correlators. Conservation of energy alone seems to suggest that such should be the case. 

Neither the mean field nor the spatial homogeneities in the fiuctuations decay away 
completely. At late times the dynamics consists of propagating wave packets in both the 
mean field and fiuctuations. This behavior may have been enhanced in the present study 
by the fact that there are no massless excitations in the spectrum and therefore there are 
kinematic thresholds prohibiting further decay. This is a general characteristic of the = 1 
theory considered here. Collectively, the characteristics of the late time evolution, 
although interesting and an improvement on the homogeneous version of the approximation 
and on the classical fields, are incompatible with a state of canonical thermal equilibrium. 
This aspect of the problem will almost certainly require the inclusion of quantum scattering 
effects 



In contrast, physical situations that are dominated by the behavior of (massless) Gold- 
stone modes -such as scalar theories with large A^- will behave quite differently. The large- A^ 
limit of these models can also be studied for spatially inhomogeneous mean fields. Relative 
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to the Hartree approximation it has the added advantage of capturing the correct (mean 
field) thermodynamics |]7,14|. Unfortunately, the study of the dynamics of symmetry break- 



ing transitions is most interesting in higher spatial dimensions, where, in order to address 
spatially inhomogeneous situations, a substantial increase in computational effort, relative 
to that of the present study, will be required. 

Our results are complementary to those of Salle, Smit and Vink 0, who studied the same 
model under additional averaging over ensembles of mean fields. This ensemble averaging is 
an additional ingredient not included in the mean-field approximation to the dynamics of the 
quantum field theory and at face value is not allowed since the mean field already carries the 
meaning of the quantum ensemble average field expectation value. It may nevertheless be the 
means to an end -the quantum thermalization of the theory- at the expense of restricting 
the mean field to particular ensembles. Our results indicate that, without invoking this 
additional ensemble averaging, the mean-field dynamics does not thermalize for a large class 
of initial conditions, including the close vicinity of thermal equilibrium. 

Because of these shortcomings the dynamical Hartree approximation is best suited to 
study situations that are not dominated by fiuctuations and that are not closely related to 
a symmetry breaking second order transition. In particular, the evolution of quasi-classical 
field configurations, such as topological defects and phase interfaces in a background of self- 
consistent fiuctuations, should prove very interesting. The hydrodynamical properties of the 
model and, in particular, the evolution of the energy momentum tensor in several situations 
are also being considered. 
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